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Quantum kinetic theory IV: Intensity and amplitude fluctuations of a Bose-Einstein 
condensate at finite temperature including trap loss 
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We use the quantum kinetic theory to calculate the steady state and the fluctuations of a 
trapped Bose-Einstein condensate at a finite temperature. The system is divided in a condensate 
and a noncondensate part. A quantum-mechanical description based on the number conserving 
Bogoliubov method is used for describing the condensate part. The noncondensed particles are 
treated as a classical gas in thermal equilibrium with temperature T and chemical potential /i. We 
find a master equation for the reduced density operator of the Bose-Einstein condensate, calculate 
the steady state of the system, and investigate the effect of one-, two-, and three-particle losses on 
the condensate. Using linearized Ito equations we find expressions for the intensity fluctuations and 
the amplitude fluctuations in the condensate. A Lorentzian line shape is found for the intensity 
correlation function that is characterized by a time constant 'yj derived in the paper. For the 
amplitude correlation function we find ballistic behavior for time differences smaller than 'yj^ , and 
diffusive behavior for larger time differences. 
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I. INTRODUCTION 

In a series of papers we have developed a quantum ki- 
netic (QK) theory with application to Bose condensation 
of cold dilute gases. In the first two papers, which we 
shall refer to as QKI R] and QKII 0, we considered a 
spatially homogeneous, weakly condensed system, where 
the interaction between the atoms was assumed to be 
sufficiently weak for quasiparticle effects to be negligible. 
In QKIII [H the theory was extended to a strongly con- 
densed gas in a trapping potential under the assumption 
that the noncondensed vapor acts as a heat and particle 
reservoir for the condensate (see also Ref. Q), a situation 
which corresponds closely to present experiments of Bose 
condensation with alkali-metal vapors P-pl| . 

In the present paper (QKIV) we will study the steady 
state, amplitude, and phase fluctuations of a trapped 
Bose-Einstein condensate at finite temperature, including 
the effects of one-, two-, and three particle losses on the 
condensate. Such a study seems particularly timely in 
view of the present interest in the dynamics and measure- 
ment of the phase of the Bose condensate (for a review 
see Ref. ]12]). Until now the discussion in the literature 
has essentially focused on phase collapse or diffusion, and 
phase revivals in the zero-temperature limit, analyzing 
the dependence of collapse and revival times on the trap 
potential, the dimensionality of the gas, atom number 
fiuctuations, and the coherent dynamics of the conden- 
sate [Ol-fsi. In contrast, in the present work we will 
study in detail fluctuations as a result of interaction of 
the condensate with a (reservoir of) uncondensed atoms. 

We will almost exclusively consider a grand canonical 
particle reservoir in this work. This particle reservoir will 
be assumed to have a constant chemical potential and not 
to be influenced by the mean field of the condensate. The 
results will therefore only be valid in the case of small 
condensate and large thermal particle numbers. In all 
other cases the simple model used here for the thermal 



particles will have to be replaced by a more sophisti- 
cated one. A more detailed discussion on the expected 
correction to the results due to particle conservation and 
mean-field effects is given in Appendix |A|. 

The starting point of the theory is to decompose the 
field operator describing the iV-particle system into con- 
densate and noncondensate parts. The division is based 
on a splitting into two energy regions where the high 
energy band is supposed to contain particles in a ther- 
mal state corresponding to a given temperature T, and 
chemical potential /i. The condensate band contains the 
actual condensate as well as quasi-particle excitations. 
A quantum-mechanical description based on the num- 
ber conserving Bogoliubov method is used for describ- 
ing the condensate part [Q. To facilitate the analy- 
sis we drop the quasi-particles and describe the conden- 
sate by a single mode with destruction operator B, and 
the spatial wave function ^Ar(x), corresponding to the 
solution of the Gross-Pitaevskii equation for N conden- 
sate particles. Elimination of the noncondensate part 
leaves us with a master equation for the (reduced) con- 
densate density operator. The physics contained in this 
equation is quite rich. The master equation accounts 
for loss/gain of particles to/from the thermal band, and 
phase-destroying but number-conserving collisions be- 
tween condensate and noncondensate particles as well as 
linear and nonlinear trap loss. 

A diagrammatic illustration of the processes described 
by the master equation is given in Fig. 1 . We realize that 
there are two types of processes. On the one hand, there 
are those which involve particles from both the conden- 
sate and the noncondensed band. They comprise pro- 
cesses which lead to a redistribution of particles between 
the two bands at rates W^{N) (condensate gain) and 
W~ (N) (condensate loss) , as well as number conserving 
scattering events of thermal particles off the condensate. 
The latter occur with a rate Roo{N) and will give rise 
to fluctuations in the condensate phase. Explicit expres- 



sions for W^ and i?oo will be given below. On the other 
hand there are several loss mechanisms at work which 
deplete the condensate [ p7|jlq ] . There is one-particle loss 
due to collisions with background gas atoms with associ- 
ated rate 71. Two particles can be lost with rate 72 (iV) 
from the condensate if two condensed particles collide, 
and one of them changes its internal state. This particle 
no longer sees the trap, and escapes. Its partner is im- 
parted with the energy difference set free by the collision, 
and is also lost from the trap. Finally, three-particle loss 
can occur with rate 73 (TV) if in a three-particle collision 
a dimer is formed. The binding energy is imparted to 
the third particle, and all of them escape from the trap. 
Note that the description of the noncondensate particles 
in terms of a thermal reservoir results in a finite occupa- 
tion of the condensate mode even in the presence of loss 
channels. 
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Fig. 1. Interpretation of the processes described by the mas- 
ter equation: (a) W'^{N) and W~{N) are the feeding and 
depletion rates of the condensate from and to the noncon- 
densed thermal band of particles, respectively, (b) _Roo is the 
rate of thermal particles bouncing off the condensate without 
changing its occupation number. 71 (c), 72 (A'^) (d) and 73 (A'^) 
(e) are the rates of one,- two,- and three-particle losses, re- 
spectively. We use a star to indicate a change in the internal 
state of a particle. The barbell represents two atoms having 
formed a molecule. 



II. MAIN RESULTS 

In this section we will give a short overview of the main 
results to be derived in detail in the remainder of this pa- 
per. As a starting point for our analysis of the particle 
number and phase fluctuations of a Bose-Einstein con- 
densate at a finite temperature, we adopt the theoretical 
description developed in the precursor papers QKI-III 



A. Fluctuation analysis 

In the limit of large condensate particle numbers, we 
may approximate the master equation [Eq. (|l^)] by an 
equation that is of Lindblad type. This has the advan- 
tage that standard techniques developed in quantum op- 
tics for the description of fluctuation properties become 
applicable [|l9[ . We have thus derived quantum stochastic 
differential equations for the condensate particle number 
N = B'^ B and the Glogower-Susskind phase operator e*"^, 
which is known to characterize phase fluctuations well in 
the limit of large occupation numbers pQ]. Lineariza- 
tion of these equations is permissible in the very same 
limit of large average occupation N = (A^), and allows 
us to work out the two-time correlation functions of oc- 
cupation number and phase. The spectra of condensate 
occupation number and amplitude fluctuations are then 
immediately obtainable by Fourier transformation. 



1. Condensate particle number fluctuations 

For the correlation function of the particle number fluc- 
tuations pi[ we obtain the following result which holds in 
the stationary limit, i.e., for times satisfying i + s 3> 77 

(1) 
with (a, 6) = {ah) - {a){h). In Eq. (Ill) 



{N{t),N{s)) = Ae-Td*-^l = cr^e-'^^l*-^!, 
27/ 



/2/2 = N{W+ + VF" + 71) + AN'^^2 + 9^^73, (2) 

with a bar denoting evaluation at A^ = iV. Note that a 
appearing in Eq. (nl) is a measure of the width of the 
particle number distribution. The characteristic time 
constant 7/ with which the particle number fluctuations 
regress is given by 



7/ = 2Ndi^W- - (72 + 373^)87V/5. 



(3) 



dfi denotes the derivative with respect to N and evalu- 
ation at N = N . The exact size of this rate depends on 
the specific experimental setup. A convenient quantity 
to assess the amount of fluctuations present is the well 
known Mandel Q parameter defined as 



Q^]iui{N{t),N{t))/{N{t))-l. 



(4) 



A coherent state would correspond to a value of Q = 0, 
while a number state yields the minimum value Q = —1. 
Assuming two and three particle losses to be insignificant 
we find for the Mandel Q parameter. 
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2/^ AT 
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(5) 



To arrive at this result the approximate expressions of 
the rates W~^ and W~ as given in Eqs. ( [l9|) in terms of 
the scattering length a, the reservoir temperature T and 
chemical potential /i have been u sed. 

The results (for details, see Sec. IV C ) can now be sum- 
marized as follows. The variance of the occupation of the 
condensate mode is proportional to kT/ ^j^, with /i^ the 
chemical potential for the condensate mode with average 
occupation N. In the Thomas-Fermi approximation ^^ 
is given by 
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(6) 
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The constants used in Eq. (pi) are defined in Sec 
The characteristic rate at which the correlation function 
drops off is roughly given by 7/ ~ 27V9jyW^~i with W 
the loss rate to the noncondensed band. 



These values have to be understood as order of magni- 
tude estimates for current experiments. Note that we 
used a value of a = 2.6nm for the scattering letngth of 
rubidium whereas recent experiments at JILA showed 
that a = 5.1nm. 



III. MODEL 



In this section we briefly describe the basic concepts of 
the quantum kinetic theory developed in Ref . ||l|-|^, p^ , p3| . 
We do not give a detailed derivation of the master equa- 
tion used throughout the paper since this can be found 
in Ref. 01. 



A. Description of the system 



2. Amplitude and phase fluctuations 

In the limit of large and well-defined average occupa- 
tion number TV of the condensatemmode, the amplitude 
correlation function {B^{t)B{s)) is well suited to assess 
phase properties of the condensate mode p^ . In particu- 
lar, the spectrum of phase fluctuations is identical to the 
spectrum of amplitude fluctuations. For the amplitude 
correlation function we obtain for t > s 



{B^t)B{s)) = 7Ve(^T^^^^™)(*"^^"''(''^(*"'^+^" 



~,,(t-s)^ 



(7) 



where 77 = (adj^fij^/jih) . As we will see in Sec. II A 3 
i?oo is negligible in the above correlation function 
[Eq. ^] for most of the current experiments. The struc- 
ture of the correlation function indicates that there are 
two distinct time regimes: 

• For 7/|i — s| <C 1, the term proportional to rj in the 
exponent is proportional to (t — s)^. This is called 
the ballistic regime. 

• For \t — s\jj > 1, the phase behaves like that of a 
system undergoing phase diffusion. A characteris- 
tic of such behavior is a linear dependence of the 
exponent on |t — s|. Note that for large time differ- 
ences we observe the legacy of the ballistic regime 
in the form of a rescaling of N to Nao = Ne^. 



3. Numerical values 

Using data from the experiments at JILA ||2^ and 
an average occupation number of iV = 25000 Rubid- 
ium atoms at a temperature T = Q.b^K in a trap with 
fx = fy — fz/VSi = 47Hz we obtain for the rates 
7/ « 2Hz, ^00 « 4mHz, M^± « 50Hz, and 77 w 800. 



The Hamiltonian of a weakly interacting Bose gas con- 
fined by a potential Vr (x) in second quantization is writ- 
ten as 



i../.3.v.t(x)(_|-V^ 



+ 1/t(x) VXx) 



+ i /"d^x /"d-''x'V'^(x)V'^(x')u(x-x')^(x')^(x). (8) 

'0(x) is the standard bosonic field operator. The two- 
body interaction potential u(x — x') is a short-range po- 
tential of the form u5(x — x') where u — A'Kfia/m with a 
the scattering length. We assume the trapping potential 



to be of the form Vt(x) = m{Lo^x 



-^ly^ 



'-)/2. 



As was shown in Ref. [p| we can obtain a master equa- 
tion for the density operator of the condensate mode by 
dividing the Bose gas into two energy regions called the 
condensate band Re and the noncondensate band i?NC- 
The boundary Eji between these two regions is chosen 
according to Ref. Q such that _Rnc is not significantly 
affected by the mean field of the condensate. The field 
operator is then written as V'(x) = '0nc(x) -I- </)(x) de- 
scribing particles in the noncondensate band and in the 
condensate band, respectively. We will treat the parti- 
cles in i?NC as classical thermalized particles character- 
ized by a temperature T and a chemical potential /z. The 
particles in Re are affected by the mean field, and they 
have to be treated quantum-mechanically. Re contains 
all the trap levels that are significantly modified by the 
mean field of the condensate. As in Ref. [g3| we will use 
the simplest possible description of the condensate band 
by assuming that only one mode, namely, the conden- 
sate mode itself, is important and all the other modes 
are negligible. The master equation we will use for our 
calculations is an equation for the reduced density oper- 
ator pc of the condensate band Re interacting with the 
bath of particles in it!NC- 



1. Condensate band Re 



2. Noncondensate band Rnc 



We use the number-conserving Bogoliubov method de- 
rived in Ref. ||l^ to describe the particles in Re- In this 
formulation we can write down the condensate band field 
operator (/)(x) in the form 



(x) = i?(iv) e^(x) + ^^i^iMzi 



Mn9r^ 



(x) 



N 



(9) 



The annihilation operator B{N — 1) brings the system 
Re from the ground state with N atoms to the ground 
state with iV — 1 atoms. The action of the operator 
B{N — 1) to the condensate is depicted in Fig. 2 . \N)p4 
denotes the ground state with TV particles in the con- 
densate. Applying the operator B to this state yields 
B\N)n = ^/N\N - 1)n-i- Note that the operator B 
changes the chemical potential of the condensate from 

MW to piAT-l- 




|W-l>»-i 



Fig. 2. Condensate ground state with A'^ particles \N)m and 
with A^ — 1 particles |A^ — 1)m~i- These two states are con- 
nected by the operator B, as described in the text. 

As shown in Ref. ||l^ Ca^(x) is the condensate wave 
function, and satisfies the Gross-Pitaevskii equation 



7^'^'^S,N + VTS,N + Nu\£,N?iN^llN^N. (10) 

2m ' ' 



In all our calculations we will use the Thomas-Fermi ap- 
proximation for the chemical potential [Eq. (p)] and for 
the condensate wave function 



C^(x) 



MAT - ^t(x) 



Nu 



for /ijv > Vt{^), (11) 



and zero elsewhere. The energy of the condensate is given 
by 



^n = 



hN^i 



N 
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(12) 



The amplitudes fm{^) and 5m (x) describe creation and 
destruction of quasiparticles. They are defined in Ref. [g| , 
but we will not need them in this paper. We will assume 
there is always a fairly large mean number of particles N 
in the condensate so that we can neglect the influence of 
the quasiparticles on ^(x) M. 



We treat the noncondensate band as a thermal bath of 
particles in thermal equilibrium. According to Ref. |3J we 
only need the phase-space density of the noncondensed 
particles i^(K,x) to calculate all the rates appearing in 
the master equation for pc [Eq. (|l7|)]. In our calculations 
we will use the classical approximation 



F(K,x) 



Jp- 



-VT(^))/kT 



(13) 



We expect corrections to the rates W^ and i?oo from 
using a more detailed description of the noncondensed 
particles. However, all the calculations presented here 
will remain valid, since they are mostly independent of 
the functional form of the rates W^ and i?oo- We only re- 
quire that it is permissible to linearize the rates W^{N) 
and Rqq^N) around the mean number of particles in the 
condensate N. In the present work we will mainly present 
results obtained by using Eqs. (|l^) and ( pO[ ) for the rates 
W^{N) and Rqo{N) which were obtained by the assump- 
tion of a grand canonical classical bath of particles not 
influenced by mean-field effects. Corrections to the sta- 
tionary state due to conservation of the total number 
of particles as well as mean-field effects and quantum 
statistics effects are estimated in Appendix ^. A fur- 
ther, more detailed, discussion of different models for the 
uncondensed bath (i.e., the noncondensate band) lies be- 
yond the scope of this paper, and will be presented in 
other work. However, we want to point out that these 
models might include canonical baths with a constant 
particle number |23], evaporatively cooled baths fP6| as 



well as baths that are continuously fed with particles p7 



3. Trap loss 

There are several processes leading to losses of con- 
densate particles from the trap. We will consider one,- 
two,- and three-particle loss with loss rates 71, 72 (A^) and 
73 (A^), respectively. One-particle loss might be caused 
by background gas particles hitting condensate particles 
or by coupling condensate particles out of the trap as 
described in Ref. |^l[. Inelastic two-particle collisions 
changing the internal properties of the particles lead to 
two particle loss. In most of the current experiments 
the two particle loss is negligible compared to the three 
particle loss caused by the inelastic collision of three par- 
ticles. In the Thomas-Fermi approximation we obtain for 
the loss rates 



71 



fJd^.\^N{.f = f, 



(14) 



72 (A^) = 



Ko 



d^x\^Nix)r^ 



K2l6-TTp.^^^V2 



ex A" 



'3/5 



(15) 



73(iV) = ^ 






(16) 



For Rubidium the constants K^ have been measured in 
|jl7| . The rates Ki have been calculated analytically in 
Refs. [|8|-|30 . 



B. Master equation 

We simplify Eq. (50) in Ref. |3| for the case that the 
condensate band consists only of the condensate mode 
alone. In contrast to Ref. [|| we keep terms including 
Rqo, and add the loss terms to the master equation (see 
Fig. 1 ). We thus obtain the following equation for the 
reduced density operator of the condensate band pc'. 



Pc = ~-[Ho,pc\ 

+ 2B^{W+iN)p,}B - [BB\{W+iN)pc}]+ + 2B{W^-(7V)pjBt - [B^B, {W^-(7V)pJ] + 
+ 2BB^{RoQiN)p,}BB^ - [BB^ BB\{RooiN)p,}] + 

+ 2S{7ipjBt _ [B^B, {7iPc}]+ + 2BB{^2{N)p,}B^B^ - Ib^B^BB, {72(^)Pc} 

L J + 

+ 2BBB{-i3{N)pc}B^ B^ B^ - \b^ B^ B^ BBB,{-tz{N)p^}\ , 



(17a) 

(17b) 
(17c) 
(17d) 

(17e) 



where Npc = \B' B, pc] /2. An intuitive interpreta- 
tion of the processes described by the master equation 
[Eq. (p^)] is given in Fig. 1 . The free evolution of the 
condensate is described by the Hamiltonian Hq, 

Ho = Eo{N) where En{N + 1) - Eo{N) = p^, (18) 

and /ijv is the chemical potential obtained from the 
Gross-Pitaevskii equation. 

Growth and depletion of the condensate due to in- 
teraction with the noncondensed particles is given in 
Eq. (|l7b|) . The growth and the depletion rates W^{N) 
of the condensate were already calculated in Ref. S us- 
ing the Maxwell-Boltzmann approximation for the phase- 
space density of the noncondensed particles. They are 
given by 



for LOz > ujx — t-^y, and 



RooiN) 



AkTp^ 



N 



9TT^h^LuluJ,N^ 



pt^/kT, 



5in(y^ 



1 — 



(20b) 



for Uz < ujx = ujy. A detailed derivation of Rqq is given 
in Appendix H. 

Trap loss is accoun ted f or by the last two lines of the 
master equation Eqs. ( |l7cj ) and (17e). Note that the only 
difference between the process including the rate W~ (N) 
and the one-particle loss rate 71 is the dependence of the 
two rates on the properties p and T of the noncondcnsate. 



W+{N) 






kT |^^^(^)| ^ (^Qa) IV. SOLUTIONS OF THE MASTER EQUATION 



w- 



-{N) = e'■^"^'^'^/'''^W+{N). (19b) 

The rate Rqo defined in Eq. (141) of Ref. [|| can be 
understood as describing the process of thermal parti- 
cles bouncing off the condensate. This process does not 
change the number of particles in the condensate, but 
does cause phase fluctuations. Using the above expres- 
sion [Eq. (p^j for F(K,x) and the Thomas- Fermi ap- 
proximation for the condensate wave function, we find 



In this section we investigate solutions of Eq. (|l7|). We 
find the stationary solution and derive a differential equa- 
tion for the mean number of particles in the condensate 
N . Using linearized Ito equations, we obtain the corre- 
lation functions {N{t),N{s)) and (B^(i)B(s)). 

We define 



55(, ■.^N{N~iy\pc\N)N. 



(21) 
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(20a) 



From Eq. (y_7|), we derive the evolution equation of the 
matrix elements, expand the terms in this equation for 
large A'^ in a Kramers-Moyal type expansion to order 
1/A', and obtain 



<?]^ = 2^{N-,.)N{W+{N - l)g5^_i - W-{N)g'^^} 



+ 2^{N-i^ + l)iN+l){W-iN + l).gXr+i - T^+ Wff^} 

+ 2v/(A^-^ + l)(iV+ 1)7135^+1 - (27V - ^)7i5]^ 

+ 2^iN -y + 1){N -v + 2){N + 1){N + 2)72(iV)g^+2 " ((^ " ^){N - v - I) + N{N - l))-f2{N)g''^ 

+ 2^{N -v + l){N-v + 2){N -v + 3){N + 1){N + 2){N + 3)73(^)3^^+3 

- ((iV - ^)(iV - ^ - l)(iV - z/ - 2) + iV(iV - 1)(7V - 2))73(7V).gJ^ - (I/t]^)^?;^. 



(22) 



The time constant 



IV/lN 



^N = {^ [W+{N) + W-iN)] + ly^RooiN) - ^ 



(23) 



determines the time scale on which the non-diagonal ma- 
trix elements decay. We have assumed that v <^ N ^ and 
therefore approximated 



{Eo{N) - Eo{N ~ y)) 



IV^N 



(24) 



A. Stationary solution 

First we neglect trap loss and solve Eq. (E3) for the 
stationary case. All g'^ with v ^ Q drop off exponen- 
tially in time, and we are left with the diagonal terms 
only. We use the detailed balance condition to calculate 
the stationary solution g%: 



gjV+l 

-9% 



and therefore 



W+{N) 
W-{N + 1) 



dli-liN)/kT 



9n oc 



N 
n=0 



,(/x-p„)/fcT 



(25) 



(26) 



Using the Thomas-Fermi approximation for the chemi- 
cal potential of the condensate, and replacing the sum 
occurring in the exponent by an integral, we obtain 



5^ oc e^(^-5A'«/7)/feT^ 



(27) 



This particle distribution has one maximum determined 
by the condition /iat — ^l, as expected from thermody- 
namics. The position of the maximum differs from the 
mean number of particles in the condensate by an amount 
of order \/N . 



1. Linearized solution 

In the case of iV 3> 1 we may linearize the solution 
Eq. (G^) around the mean number of particles in the con- 
densate TV. We approximate the distribution in Eq. (127 
by a Gaussian 



-9% 



2iTa 



(28) 



The width a of this Gaussian is given by 



IbkTN 



2fj. 



(29) 



N 



In the Thomas-Fermi approximation a scales with the 
mean number of condensate particles like 7V^/^°. The 
difference between the linearized and nonlinearized solu- 
tions can therefore only be seen for very small conden- 
sates. Figure 3 shows a comparison between the Gaussian 
approximation and numerical solutions obtained from 
Eqs. (|l^) and (|3^). As expected both solutions agree 
very well with each other even for a mean occupation of 
the condensate of only about N « 500. Note also that 
the same result for the variance may be obtained from 
statistical mechanics. Thus the restriction fij^/kT > 1 is 
not necessary for this result to be valid. 



9n 




Fig. 3. Stationary particle distribution in the condensate. 
The trap loss is assumed to be zero. fi/kT — 0.05 and 
^N=i/kT — 0.004. The solid line represents the numerical 
stationary solution of the Fokker-PIanck equation [Eq. (p0[)] 
and the detailed balance solution of the master equation 
[Eq. (h_7|)]. The dashed line depicts the Gaussian approx- 
imation [Eq. (ESI)]. The trap frequencies are chosen to be 
fx~fy = fz/yS — 47Hz. The calculation is done for Rubid- 



2. Inclusion of trap loss 

In our model the bath of thermal atoms is not depleted 
by the interaction with the condensate. Experimentally 
this can be achieved by replenishing the reservoir by some 
mechanism, or by doing the experiment so quickly that 
the number of particles lost from the reservoir can be 
neglected. Furthermore, the calculations presented here 
for constant T and /z remain valid as long as the heat 
bath parameters fi and T change slowly compared to the 
time scale of the condensate dynamics. For the diago- 
nal matrix elements g'^, we therefore obtain a stationary 
solution different from zero even if we include trap loss. 
Keeping only the leading order terms in N of Eq. (22), 
we immediately find the Fokker-Planck equation for 5^ 
to be 



9l 
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{fiiN)9%} 
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dN' 
where we have defined 



2dN'^ 



{f2{N)g%} (30) 



and 



f^{N) = -2NW+{N) + 2N{W-{N) + 71) 

+m^j2{N) + 6N^j3{N), (31) 

/2(iV) = 2NW+{N) + 2N{W-{N) + 71) 

+8N^j2{N) + 18N^-f3{N). (32) 



The Fokker-Planck equation is valid as long as iV ^ 1 
and a :s> I. 

We approximate the solution of this Fokker-Planck 
equation by a Gaussian, and obtain the following equa- 
tion for the mean number of condensate particles iVioss: 



/i - 0. 



(33) 



The width of the Gaussian ctioss is approximately given 
by 



CTlo 



/2 



29^0./! 



(34) 



Using the assumption fij^^^^JkT <C 1 Eq. ( p3| ) can be 
solved analytically by approximating 

W+{N) « W+ = 4"^(°fc^)' e2MAT^ (35) 



W-{N) « W-(N) = W^{1 + (^ + ^)e-^'A^. 

(36) 

For the mean number of particles in the condensate, we 
obtain 



A^loss = 



'373(1) + M?W^a+e-A'/fe^/2 V ^73(1) + ^JilW+e-^^/'''r/2) 373(1) + /ifW+e-^A^/2 



5/2 



(37) 



Here /xi = fj.]\[^i/kT. Trap loss decreases the number of 
particles in the condensate. Also, the width of the parti- 
cle distribution is decreased by the nonlinear loss. Both 
of these effects are well known in nonlinear optics [Q . 
Figure 4 shows the effect of trap loss on the mean number 
of particles in the condensate for a given T and /i for the 
parameters measured in Ref. p7| . The dominant con- 
tribution to the trap loss comes from three-particle loss, 
while one,- and two-particle losses are almost negligible. 
In the following sections we will omit the subscript loss 
since all the calculations will be done for finite trap loss. 



cles in the condensate for rubidium. The trap frequencies are 
fx = fy ~ /z/\/8 — 47Hz. The temperature of the thermal 
cloud is chosen T — 0.5^iK. One,- two,- and three-particle 
losses are given by the dotted, dashed, and double dotted 
lines, respectively. The solid line accounts for all three kinds 
of loss. We have used Ki = l/70s"\ K2 = 10"^^m^/s and 



B. Nonstationary solutions 
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Fig. 4. Influence of trap loss on the mean number of parti- 



From Eq. (O) we find the evolution equation 
the mean number of condensed particles TV. 

N = 2W+iN)iN + 1)- 2W~{N)N 

-2jiN - 472(7V)iV2 - 673(^)^^ 



m for 



(38) 



By comparing Eq. ( pq ) with the results obtained in 
Ref. ||l^ we find numerical values for the constants K^. 
Equation (^) was investigated in Refs. [^,Q for the 
case ji — 0. Here we only want to show that trap loss 



does not change the general behavior of the growth of 
the condensate. In Fig. 5 we show a comparison between 
the growth curve with and without trap loss. For the 
parameters chosen the time scale of condensate growth 
is influenced only slightly by trap loss. 



N 




Fig. 5. Comparison between the growth of the condensate 
with trap loss (solid curve) and without trap loss (dashed 
curve). The parameters are chosen as A'l — l/70s~^, K2 — 
10"^^mVs, and K3 = 6 ■ 10"*^m''/s |0. Parameters for the 
thermal particles are T = 0.5 jjK and /i = 3- 10~^^ J. The trap 
frequencies are chosen to be fx ~ fy = fz/VS = 47Hz. The 
calculation is done for rubidium. 



C. Correlation functions 



In this section we calculate the intensity and amplitude 
correlation functio ns for a condensate in the stationary 
state found in Sec. IV A, 



1. Ito equation 

First we will show that if we restrict ourselves to sit- 
uations where a large number of particles occupies the 
condensate and the density operator is almost diagonal 
in the basis \N)n, we can approximate Eq. ( |l7| ) by a 
master equation of standard Lindblad form. To do so, 
we consider terms of the form 



2CHF{N)p}C~[CC\{F{N)p}] 



(39) 



define the operator D = \J F{N)C, and compare Eq. (|3 
with 



2D''pD-[DD'',pl 



(40) 



for a matrix element g'^. Since we assume the trap loss to 
be a small effect compared to the interaction of the con- 
densate with the thermal particles, the terms C^CF{N) 
of the master equation ( p7| ) are of order W~^N. The only 
exception are the terms involving Rqq |34[ |. Therefore, 
we find for all the terms in our master equation that the 
difference of Eqs. (p9| ) and (^0|) for a matrix element g'^ 
is of the order of NW~^{N){i'/N)'^. We will linearize 
the equations and neglect all terms of order W'^(N)/N. 
Therefore, we approximate the master equation (pTJ ) by 



replacing all the terms of the form Eq. ( p9D by expres- 
sions of the form of Eq. ([iol). This enables us to write 
down the Ito equation straightforwardly for the evolution 
of the system operators in the Heisenberg picture. Note 
that the noise terms appearing in the Ito equations do 
not have a direct physical interpretation. We only need 
them to have a mathematical equivalence between the 
solutions of the master equation and the solutions of the 
Ito equation. The Ito stochastic equation for an operator 
X in the Heisenberg picture reads 



dX = i.^[Ho,X]- \X,Vw+BB^Vw+]+ + 2Vw+BXB^Vw+- [X, JW- + -fiB^ bJW- 
l n 



-71J 



'RnaBB^XBB 



2^W- + jiB^XBy/W- + 71 - [X, y/R^a{BB^fy 

-[X,^B^B^BB^^+ + 2^B^B^XBB^~ [X,^B^ B^ B^ BBB^]+ + 2^B^ B^ B^XBBB^^t 

- {[X, V2W+B]dCi - [X, BW2W+]dCl^ - ([X, ^2{W- + -fi)B^]dC2 - [X, B^2{W- + -n)]dCl') 



[X, y/2R^oBB^dC3 - [X, y^2R^oBB^]dCl 
([X, y/^f^B^ B^ B^]dC5 - [X, BBB./2^]da 



W- + 7i)B^]rfC2 - [X, By/2{W- + 71 J 
[X, y^B^B^jdCi - [X, BBy/2^]dCl' 



5 1 > 



(41) 



where dCi are Ito noise increments. The only expectation 
values that are different from zero are 



{dC^{t)dC\{t)) =dt. 



(42) 



Note that all the rates appearing in Eq. (41) depend on 
the number operator N and, that therefore we use rela- 
tions like [for example, for W~{N)\ 



BW-{N)^W-iN+l)B^ (w-{N) + ^^^^4^]b 

= {W- +djyW-)B (43) 

to calculate commutators between these rates and B. 

2. Intensity fluctuations 



We define the operator 



Sip 



N- (AT) 



(44) 



where we omit the time dependence whenever this can 
be done without causing confusion. For SIb we obtain 



dSIs = —jiSlBdt + dCi, 
where 
dCi = (V2W+BdCi + BW2W+dCl~ 



(45) 



-^2{W-+ji)BUC2 - B^/2{W~+-fi)dCl 
-2y/¥f2B^BHCi - 2BB^/2^dCl 



-3^/2^B^ B'' BUC5 - iBBB^/^ildCl 



1 

/n' 

(46) 



We expand around N and obtain 
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II = df,h « 2Ndj^W~ + -72 iV + —i:iN\ (47) 

which is now only a function of the expectation value of 
N . Keeping only these highest-order terms, we can solve 
Eq. ( ^ ) and obtain 

5lB{t) = 5lB{0)e-^'' + I e-^^(*-*') dCi{t'). (48) 



Using {dCi{t)dC\{t)) — f2dt/N and considering only 
t + s 3> 77 we obtain for the correlation function 



(iV(i),7V(s)) = Ae-7.l*-l 
27/ 



cr^e-^^l*-"!. (49) 



Note that the operator SIb is of order 1. For 72 = 73 = 
we find for the Mandel Q parameter 



5kT 



1. 



(50) 



This means that as long as there is a significant ther- 
mal bath, Q will always be larger than since in this 
case fj,f^/kT < 1 holds. However, since the result for 
Q also follows from statistical mechanics the restriction 
fij^/kT < 1 is not necessary for this result to be valid, and 
we obtain sub-Poissonian statistics for ^f^/kT > 5/2. 



3. Amplitude fluctuations 

We use the Ito equation introduced above to calculate 
the phase fluctuations. In particular we calculate the 
correlation function 

{BHt)Bis)). (51) 

First we simplify Eq. (BTl) ior X = B and find 



dB = -jBBdt 
where to leading order in N 



IB 



h 



25^°" 



NSIf 



dCp 



idj^jij^ 



(52) 



II 32 J?oo \ 
2N 125N ) 



(53) 



The terms of order W^ /N have been neglected. Expres- 
sions like Eq. ( p3[) appear in optics in connection with 
the Kerr effect [B5|. The noise term reads 



dCE 



f Vrn^ 



OmW^ 



V2W+ 



B^B ] dC\ 



jn\ 



V2W+ 



BBdCi + -yj2{W- +-11) 



dNW- 



V2(M^--f7i) 
+ ^BBBdCl 



-.BBdCl 



dNRoo 



B^B - V2-R00 ) B (dC3 - dCl 



'W2 



/2^ 



y2Ro{) 
B^B^B^B - 3V^B^B^ ) dCr^ 



OnW- 



^2{W-+^i) 



-.BB^ dC2 



/2^ 



B^B^S - 2^/2^2B^ ) dC, 



/272 



BBBBdCl 



(54) 



We define the operator $ = 
lowing equation for $ 



, ^ B and find the fol- 

V-B + -B + 1 



d$ = -jBdt + dC^ T=d6lB * (55) 

V 2VN J 

Keeping only the leading terms, we obtain for dC^ 



and find 



V2Ro 



B^B - V^R^) {dC3 - dC|) (56) 



dC4t)dc^{t) 



16 
25' 



Rn 



NSlB^^^]dt. (57) 
125N ' ' ' 



Equation ( p5| ) can be treated as a c-nuniber equation for 
$, since i times the noise term in this equation has the 
properties of a classical noise term. We define an opera- 
tor $ = e**^ and obtain, for its increment, 



' h 



dj^Uff 



NSIb dt - idC^, (58) 



In this equation the intensity noise dCj is independent 
of the phase noise dC^, so that the equation for cj) is very 
similar to the equations appearing in the phase diffusion 
model for a laser with colored noise ||3^ . We may use the 
same method as used in quantum optics to calculate the 
correlation function {B^{t)B{s)), and obtain 



(B^{t)B{s)) = 7Ve('^"*t^™)(*-"^ 
^ (-(^%^)'(7.(*-)+e-^<-=)-i)) 



(59) 



for i > s in the stationary state. In the range of validity 
of our approximations the result for this correlation func- 
tion Eq. (p9) agrees with the result obtained in QKIII 
[Eq. (184) of Ref. Q]. If we were to keep the terms of 
0(l/ViV) in the noise terms the phase noise would no 
longer be independent of the intensity noise and the so- 
lution to Eq. (Isl) would not be so easy to find. The cor- 
relation function ( p9[ ) depends only on the time difference 
t — s. The spectrum is found by a Fourier transformation 
in the time difference t — s 



Siu) 



1 



dAte-''^'^^B^t + At)B{t)). (60) 




The mean number of particles therefore is iV = 23800. Plot 
(b) shows the spectrum for T — O.Q/iK and fi = 2-10"''^ J, and 
therefore a mean number of condensate particles N — 8640. 
Fig. 6 shows a comparison between the analytic result 
[Eq. (|59|)] and the direct numerical solution of Eq. ^\J\). 
For a large number of condensate particles N, the results 
agree very well with each other. In case of small particle 
numbers N the linearization used to obtain the analytic 
result shifts the curve compared to the numerical result. 
Even so, the shape of the solution is very well approx- 
imated by the analytic formula. The spectrum S{lo) is 
expected to be of Lorentzian shape around its maximum 
value. Further away from the center, the shape becomes 
Gaussian. However, for the parameters chosen in Fig. 6 
the Gaussian part dominates. 



V. CONCLUSIONS 

In this paper we have calculated the correlation func- 
tions for amplitude (phase) and intensity fluctuations of a 
Bose condensate due to interactions with a heat and par- 
ticle reservoir, representing uncondensed atoms at finite 
temperature. The present analysis is valid for a strongly 
condensed system confined in a trapping potential, ig- 
noring contributions from quasiparticle excitations [pj. 
Finally, we point out that the present theory is readily 
adapted to a class of highly interesting problems, such 
as the study of decoherence in Josephson-like situations, 
where two trapped condensates are brought into contact 
and the quantum dynamics of the relative phase of the 
two condensates is observed [O 37 . 
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Fig. 6. Spectrum of the amplitude fluctuations S{uj) as de- 
fined in Eq. (pOl) against uj for rubidium. The trap frequencies 
are chosen to be fx = fy = fz/V^ = 47Hz. The numerical 
solutions are given by the solid lines, and the dashed lines rep- 
resent the analytical results. In (a) the parameters of the ther- 
mal cloud are chosen to be T = 0.25fiK and p = 3 • lO"^'^ J. 



APPENDIX A: STATIONARY SOLUTION FOR A 
CANONICAL BATH OF PARTICLES 

We want to investigate the fluctuations in the num- 
ber of particles of a Bose-Einstein condensate, assum- 
ing the system to be in the canonical ensemble. We 
will include interactions between the Bose particles in 
our analysis, and investigate their effects on the conden- 
sate fluctuations. On first glance one might expect that 
it would be satisfactory to account for the interaction 
effects by just including the mean field of the conden- 
sate (in the Thomas- Fermi approximation). However, in 
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this approach fluctuations in the size of the condensate 
lead to an unrealistically large shift of the energy levels, 
which prevents any condensation pq |. The inclusion of 
the mean field of the thermal density of particles reduces 
the shift of the energy levels. As long as fluctuations 
in the chemical potentials of the thermal cloud fi and 
the condensate ^n are small compared to the energy gap 
between the condensate energy and the first excitation 
energy, fluctuations will not lead to degeneracy in the 
thermal cloud. The excitation energies therefore depend 
mainly on the total number of thermal particles M and 
the number of particles in the condensate N. We will 
assume that the eigenenergies of the excited levels de- 
pend on M and N but are independent of n„i, i.e., the 
microstate of the system. As shown in Ref. |3^, the 
density of particlelike states is much larger than the den- 
sity of quasiparticles. Thermodynamic quantities of the 
Bose gas are therefore mainly determined by the par- 
ticlelike states, which allows us to use the Hartree-Fock 
approximation for describing the interactions in the ther- 
mal cloud. 



1. Excited modes 

We denote the energies and the wave functions of the 
excited states by e^ and ^m(x), respectively. The oc- 
cupation of these levels is rim- In the Hartree-Fock ap- 
proximation 13911 , the effective potential for the thermal 
particles Kiff (x) is given by 

Kff(x) = Vri^) + 2uN |^jv(x)|' + 2un(x), (Al) 

with n(x) the density of the noncondensed particles. 



2. Weakly interacting Bose gas in the canonical 
ensemble 

The density operator of a Bose gas in the canonical 
ensemble is given by 



1 



-f}H 



(A2) 



where /3 = 1/fcT, with T the temperature of the sys- 
tem. The partition function Zc is given by Zc = 
tr{e~''^}. We want to investigate the properties of a 
Bose condensate in the canonical ensemble. The eigen- 
states {^Ar(x), i^,„(x)} and the corresponding eigenener- 
gies {Eq, em} depend on the number of condensate parti- 
cles N and on the number of particles out of the conden- 
sate M = J2m '^™ ■ '^^^ total number of particles A^tot is 
constant, 



Nt, 



N + M ^ const. 



(A3) 



The state of the system with N particles in the con- 
densate and Hm particles in the levels e,„ is denoted 



by |n) where n = {N,ni, ...,nm, ■■■}■ We can there- 
fore write for the matrix elements of the density operator 
p(n) = (n|pc|n) 



p(n) ex e 



-l3EoiN,M)-l3j2 >^m{N,M)n„ 



(A4) 



As can be seen from Eq. (A4), the condensate energy 
Eo{N, M) and the excitation energies €m{N, M) are func- 
tions of the number of condensate particles A'^ and the 
number of noncondensed particles M. In our calcula- 
tions we will assume that the influence of the number 
of noncondensed particles M on the condensate energy 
is negligible, since the number of noncondensed particles 
that are inside the condensate region is much smaller 
than the number of condensed particles. Moreover, the 
width of the condensate particle distribution is only in- 
fluenced by the change of the number of noncondensed 
particles in the condensate region due to fluctuations in 
M. The other interaction effects i.e., (i) the influence 
of the condensate mean field on the excited levels, (ii) 
the influence of the mean field of the thermal cloud on 
the excited levels, and (iii) the influence of the conden- 
sate mean field on the energy of the condensate, will be 
included in our calculations. 



3. Particle number distribution of the condensate 



Since we are only interested in the probability of find- 
ing A^ particles in the condensate we want to find 



p{N,M)= Y,PiN,nm) 



(A5) 



summed under the constrained M ~ J2m "■™- ^^ '^^^ ^^ 
this summation by introducing a contour integral writing 



p(A^,Af)oce-'3^''(^^^'^) 



(A6) 



1 f dz 



-M 



27ri ,/c ^ 



HE 



■^ 11 2. « 

m n™ —0 



~/3e„(Ar,J\/)n„ n„ 



and integrating using the method of steepest descents to 
obtain 

p(Ar,M)cxe-'^^"(^^*^)-(^'^+i^'"("")-S™'"(i-""'^"''""*"'"') 

\ -1/2 

^ • (A7) 

z depends on Ai and A^, and is given by the solution of 

1 




m 

By defining 



,/3e„(Af,Af)2-l _ 1' 



(A8) 
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M 
F{N, M) = E^{N, M) + — ln(z) 



-ij^lnfl-z-e 



-Pe„^{N,M) 



C{N, M) 



ln(9A/ln(z))-ln(z^ 



(A9) 

l3dMern{.N,M) 
p/3e,„(Ar,M) _ £ 



+ z- 



T 



l^M^^il^^^AU (Aio) 



(e'3'-W*^)-z)^ 



and using Eq. (A8), we may rewrite Eq. (A7) 

p{N, M) ex e-'3^(^>^)+cW*^). (All) 

C{N, M) is a small logarithmic correction to /3F{N, M) 
that we neglect in our further calculations. Consistent 
with the neglect of C{N,M), we approximate M + 1 by 
M in Eq. (|A8|) and use the equations 



m=j:. 



1 



./3£„(Af,M)^-l _ 1' 



and 



p(iV,M)(xe-^^(^'^^). 



(A12) 



(A13) 



F{N, M) is the Helmholtz free energy of the system. 
The chemical potentials of the condensate and the ther- 
mal bath are therefore given by the partial derivatives of 
F{N, M) with respect to N and M, respectively. 



a. Stationary solution 

The chemical potential of the condensate in the canon- 
ical ensemble fidN, M) is given by 



^i,{N,M)^dNF{N,M) 



= OnEo [N, M)+2_^ ^,^ _ - — , ( A14) 



and the chemical potential of the noncondensed particles 
in the canonical ensemble iiv{N, M) reads 



fi,iN,M)^dMFiN,M) 
ln(z) 



/3 



E 



dMEoiN,M) 
zdMem{N,M) 



(A15) 



If the total number of particles is fixed, the condensate 
particle distribution will reach its maximum value for 
N = N and M = M = Ntot - N, where N can be found 
by solving 



The width of the stationary canonical distribution cTcan 
is given by 

dean = {Pdr^ {t^c{N, TVtot -N)^ fi^N, iVtot - N))y'/^ . 

(A17) 



In the above equations (A14) and (A15), the terms in- 
cluding derivatives of e™ {N, M) are corrections to the 
chemical potentials due to the shift of the excited energy 
levels with changing M and N. The term dMEo{N,M) 
in Eq. ( Al5| ) is the correction to the chemical potential 
of the thermal bath due to changes in the condensate en- 
ergy. By assuming Eq{N,M) « Eo{N), we will neglect 
this term. 



b. Method of solving for p{N, M) 

For simplicity we will assume that the energy and wave 
function of the condensate are only functions of the num- 
ber of condensed particles, and use the Thomas-Fermi 
expressions ( \\/\\ j and (|l2|), respectively. The influence of 
the thermal particles on the energy and wave function 
of the condensate will become important at temperature 
close to Tc and may therefore only be neglected at very 
low temperatures [ESl . 



We replace the sum in Eq. (A12) by an integral 



E 



deg{e) 



(A18) 



where g{e) is the density of states and use the semiclas- 
sical expression 

(A19) 

to obtain, for a spherically symmetric effective potential 
and given M and N, 



M ^ Att 



dr h{r) 



(A20) 



(^^T^^'\^rr^drG.J^-^-H^ 



Utt^V 



'3/2 



kT 



Here Ga{x) = S^i '^"'^^""'^ denotes the Bose func- 
tions. By solving this equation numerically we obtain 
z as well as the density of noncondensed particles n[r). 
Given a fixed total number of particles Ntot and a tem- 
perature T of the system, we find the particle distribution 
function p(A^, Af) from Eq. ( A13 ) by replacing 



-^Infl-ze^ 



-/3e™(Af,Af) 



(A21) 



2mkTV''^ TT '■°° 



ttTi 



'drG 



5/2 



K.ff(x) 
kT 



- ln(z) 



/i,„(7V, A^tot - N) = HciN, Ntot - N). (A16) in the expression Eq. (M) for F{N, M). 
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c. Solutions 

We want to compare p{N, M) with the solution for a 
grand canonical bath of thermal particles not influenced 
by mean- field effects as given in Eq. (28). To distinguish 



between the effect of choosing a different thermodynamic 
ensemble and the effect of including the change of the 
mean field due to condensate fluctuations, we will plot a 
third distribution obtained from the canonical ensemble 
with a fixed mean effective potential T4ff (x) that is given 
by 

Feff(x) = Ft(x) + 2uN |e^(x)|' + 2Mn(x), (A22) 

where fi(x) is the thermal density obtained for N par- 
ticles in the condensate, and A^tot — N particles in the 
thermal cloud. 
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Fig. 7. Condensate particle distribution for (a) T = 0.35/iA', 
iV = 1.595 ■ 10^ and N/Ntot « 10%, and (b) T = 0.2fiK, 
N = 1.910 • 10^ and N/Ntot « 33%. Both plots are for Rb 
and a trap frequency of fx = fy = fz ~ 66Hz. The solid line 
shows the solution for the canonical ensemble including the 
mean-field efl^ects, the dashed line is the canonical solution us- 
ing the mean effective potential K:ff(x), and the dash-dotted 
line is the grand canonical result. 



In Fig. 7 , we plot the different particle distributions. 
For the three curves in each set the number of particles 
in the condensate N is the same. Since the inclusion of 
the mean-field effects shifts the chemical potential of the 
noncondensed atoms, the two canonical solutions are not 
plotted for the same total number of particles. If we were 
to plot the canonical solutions with and without inclu- 
sion of the mean-field effects at the same total number 
of particles A'tot , their maxima would be shifted against 
each other. As can be seen from Fig. 7 , the stationary 
state of a Bose condensate in the canonical ensemble is 
different from the stationary state of a Bose condensate 
interacting with a bath of particles in the grand canon- 
ical ensemble. When the ratio N/Ntot ~ 10% this dif- 
ference is not significant, but for larger condensates, i.e., 
N/Ntot > 30%, the difference might become substantial. 
As shown in this appendix, the model used to describe 
the thermal bath of atoms may have some influence on 
the particle fluctuations in the condensate. We also ex- 
pect the intensity relaxation rate 7/ to be influenced by 
the particular model used for the noncondensed particles. 
However, a further more detailed analysis of such effects 
including a treatment of the excitation modes as given in 
Ref. p^ lies beyond the scope of this paper, and will be 
presented elsewhere P5|,p6[. 



APPENDIX B: CALCULATION OF Roo{N) 



We show how the expressions (|20|) are derived from 
Eq. (141) inRcf. g. 



1. Simplifying the integral 

i?oo is defined by 



RooiN) ^ ^" , j SvL I d?^i I cPK2 /"d^k /"d3k'<5(Ki-K2-k + k')i^(Ki,u)(l + F(K2,u)) 



(27r) 
Wo(u,k)Wo(u,k')<S(Aa;i2(u)), 



where 



M^o(u,k) = ^y'd'Ve^(u+^)e^(u- ^) 



ik-v 



(1) 



(2) 



and Aci;i2(u) accounts for energy conservation. We assume that within the range of the condensate F(K, u) is 
constant, and that the factor 1 -f F{K.2^ u) w 1. Integrating over K2 and defining Q = k — k' yields 



Rm{N) 



4u^ 2m 



y'dvy'd3„^|dvy'd3^^^(u+^)c^(u-:^)G(u+^)e^(u-^)e^'^'(-+-')+'Q\ (s) 
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By integrating over k' and over v' and defining 



G\Q) = 



1 



(27r)3 



^N{n~-) 



2 



Civ(u + 



JQ-v 



(4) 



we obtain 



RooiN) = 



Au^ 2m 



(27r)5?i2 n 
^(Ki,0)G2(Q). 



I £Ki I (fiq, 5{q, ■ (Q + 2Ki)) 

(5) 



Now we approximate Q + 2Ki k, 2Ki, since G^(Q) is 
sharply peaked around Q = compared to the width of 
F(Ki, 0). Replacing the notation Ki by K, choosing the 
Kz axis in the direction of Q, we obtain 



Rno{N) 



Au^ra 
'(2nf¥ 

i'^IkI + k'^ + kI) 



(fq, f dK, I dKy f dKJiQK, 



+M)/fcTg2(Q)^ 



(6) 



Here Q denotes the modulus of Q. Then we do the inte- 
gration over K to obtain 



Roo{N) 



iu irm kT 
(27r)5?i5 



.f^/kT 



,3r.G'iQ) 



rf-^Q- 



Q 



(7) 



So far we have not assumed a particular form of the trap- 
ping potential or the condensate wave function. These 
properties are contained in the function G(Q). 



and therefore 



G(Q) - ^^^l^^x |^^(x)|^*Q-. 



(9) 



Using the Thomas-Fcrmi approximation for the conden- 
sate wave function, defining Q'^ — Qx/y/o,, Q'y — Qy/vb 
and Q'z = Qz/\/c, and changing the variables x' ~ x^/a^ 
y' — yVb and z' — z^fc^ simplifies the integral to 



G'(Q) = 



(27r)3/2Va6c7Vu 



d'yl{yiN - X 'y^ ■ 



'2n iQ'x' 



X </ijv 



(10) 



Integrating in spherical eoordinates we obtain 

2 C^/P-N 



G(Q) 



dr {r^N — r ) sm{Q'r) 



y/2'KabcNuQ' Jo 

VabcNuQ ° 

^y/JJ^Q' cos{Q' y/JI^)] , 



(11) 



where Q' = JQ'^ + Q'^ 



C. Integration over Q 



B. Calculating the function G(Q) for the harmonic 
oscillator 

The trapping potential is of the form Vt(x) — ax"^ + 
by^ + cz"^ with a = muj'^/2, b = mu)y/2 and c = muj'^/2. 
We change the variables by defining x = u + ^ and 
y = u — ^ and get 

(8) 



Last we will now find an expression for the integral 

j7=/d3Q^!M. (12) 



We make use of the fact that G(Q) only depends on 
Q' by replacing d^Q = \fahc d'^Q' using spherical co- 
ordinates, and change the variable Q' — > ^JJinQ' = T. 
The modulus Q can be expressed in terms of Q' as 
Q == Q' \Ja sin^ d cos^ </> + 6 sin^ d sin^ + c cos^ Q. We 
obtain 



^ 



8/i^ 



dB- 



sin(6l) 



TTVahcN'^u^ Jo Jo \J a sin^ cos^ + 6 sin^ sin^ (i> ^ c cos^ i 



y dT— [(3 - T^) sin(r) - 3T cos(T)] ' 



(13) 



Next we integrate over T, define x — cos 0, o? — 
c — acos^ (f) ~ bsir? (f), and 0^ — (acos^ (f> + bsvc? (/))/a^, 
and write 



0.= 



IJ-% 



dnVabcN^u^ 



dx 



1 V^^TW 



2f^% 



^■n\/ abcN'^v? Jo 



2tt 



arcsinh(/3 ^) 



(14) 



In the case a^ < we take the modulus of a^ and re- 
place arcsinh -^ arcsin. To further evaluate this inte- 
gral we now assume that a = b. Then a^ = c — a and 
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/3^ = a/{c—a). The integration over cj) then yields a f acto r 
of 2tt , and we obtain the expressions given in Eqs. (20a) 
and ( [20b| ) . Roo{N) therefore depends on the number of 
particles in the condensate as N~'^/^. 
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